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SOLUTIONS OF MULTI-COMPONENT FRACTIONAL SYMMETRIC SYSTEMS 


MOSTAFA FAZLY 


Abstract. We study the following elliptic system concerning the fractional Laplacian operator 

(—A) Si Ui = • • ,Um) in R n , 

when 0 < < 1, u* : R n —> R and Hi belongs to C 1 , 7 (R m ) for 7 > max(0,1 — 2min{sj}) for 1 < i < m. 

The above system is called symmetric when the matrix H = (djHi(ui, • • • , Wm))7j = 1 i s symmetric. The 
notion of symmetric systems seems crucial to study this system with a general nonlinearity H = (L6)2Li- 
We establish De Giorgi type results for stable and H -monotone solutions of symmetric systems in lower 
dimensions that is either n — 2 and 0<Sj<lorn = 3 and 1/2 < minf-s^} < 1. The case that n — 3 
and at least one of parameters S{ belongs to (0,1/2) remains open as well as the case n > 4. Applying a 
geometric Poincare inequality, we conclude that gradients of components of solutions are parallel in lower 
dimensions when the system is coupled. More precisely, we show that the angle between vectors Viti and 
Vuy is exactly arccos (\djHi(u)\/djHi(u)). In addition, we provide Hamiltonian identities, monotonicity 
formulae and Liouville theorems. Lastly, we apply some of our main results to a two-component nonlinear 
Schrodinger system, that is a particular case of the above system, and we prove Liouville theorems and 
monotonicity formulae. 
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Keywords: Nonlinear elliptic systems, fractional Laplacian operator, Hamiltonian identity, monotonicity formula, symmetry 
of entire solutions. 

1. Introduction 

We examine the following system of nonlocal elliptic equations for 1 < i < m, 

(1.1) (-A) Si Ui = Hi(m, ■ ■ ■ ,u m ) in R", 

when (— A) Si stands for the fractional Laplacian operator and Ui : R™ —> R. We assume that each parameter 
Si belongs to (0,1) and H = (Hi)^L 1 is a sequence of functions that each component Hi belongs to C 1 ’ 7 (R m ) 
for 7 > max(0,1 — 2min{si}). We also assume that the nonlinearity H satisfies 

(1.2) diHj(u)djHi(u) > 0 when djHi(u) := f or 1 < i < j < m. 

Ollj 

More recently, the fractional Laplacian operator (—A) Si has been of great interests in the literature and 
various properties of the operator are explored. There are different mathematical approaches to define this 
operator. The fractional Laplacian operator on R" can be defined as a pseudo-differential operator using 
the Fourier transform 

(1.3) (-To^(C) = ICI 2 s ^(C), 

when the hat operator is given by 

(1.4) w{C) = (27r) _ri / 2 f w{x)e~ l< ^' x dx. 

J R» 

Suppose that each Ui £ C ,2<T (R") for a > Sj > 0 and J R „ < oo. Then, the fractional Laplacian 

of Ui can be also defined as 

(1.5) (“A )«»,(X) := P.V. / 

J R n \X-Z\ n + 28 ' 

The author gratefully acknowledges Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery 
Grant and University of Texas at San Antonio Start-up Grant. 
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for every x G 1" where P.V. stands for the principal value. It is by now standard that the fractional Laplacian 
(—A) s operator where s is any positive, noninteger number can be denoted as the Dirichlet-to-Neumann map 
for an extension function satisfying a higher order elliptic equation in the upper half space with one extra 
spatial dimension, see [11] by Caffarelli and Silvestre. To be mathematically more precise, let .s, G (0,1) and 
ui G C' 2o '(R n ) fl L 1 (R”, (1 + \z\) n+2si dz) when a > Si for all i =!,■■■ , to. For x = ( x , y) in R” +1 , set 

r y 2,Si 

(1.6) Vi(x) = / Pi(x - z,y)ui(z) dz when Pi(x,y) = p n , Si - 


and p n ,si is a normalizing constant. Then, v = (ui)’A 1 satisfies the following system 

1 div(y 1_2si Vt>j) =0 in R” +1 

(1.7) J Vi — Ui on <9R+ +1 , 


— lim y 1 2si d y Vi = d Si {—A) Si Ui on i9R!( +1 , 

for the constant 

, r(i - Si) 

(L8) dsi 2 2 ^- 1 r( Si )’ 

and Vi G C 2 (R™ +1 ) fl C(R" +1 ) and y 1 ~ 2si d v Vi G C(R" +1 ). In short, the extension function v = (i’i)^L 1 
satisfies the system 

, f div(y ai Vvi) = 0 in R " +1 = {iG R", y > 0}, 

\ -lim^o y ai d y Vi = d Si Hi(v) in 5R ” +1 = {x G R”, y = 0} , 
when a,i = 1 — 2sj. For the sake of simplicity, we fix the following notation. 

Notation 1.1. Suppose that 0 < Si < 1 for 1 < i < m. Then, s* = mini<j< m {sj} and s* = rnaxi<j< m {si}. 
We now define the notion of stable solutions by linearizing system (1.9). 

Definition 1.1. A solution v = (vi) r T_ 1 of system (1.9) in R " +1 is called stable, if there exists a sequence 
<t> = (^i)KLi C'°°(R" +1 ) satisfying 

( di v(y ai Wc/)i) = 0 in R” +1 , 

\ ~ lim y ^. 0 y ai d y <j>i = d Si Y,™=\djHi{v)<t>j In 9R!( +1 , 
when djHi(v)(j)j<t>i > 0 for all i,j and each <j>i does not change sign for all i = 1,...,to. 


In 1978, De Giorgi [20] conjectured that bounded monotone solutions of the Allen-Cahn equation 

(1.11) Au + u-u 3 = 0 in R”, 

are one-dimensional solutions at least up to eight dimensions. This conjecture is known to be true for 
dimensions n = 2,3 by Ghoussoub-Gui in [36] and Ambrosio-Cabre in [2], respectively. There is an example 
by del Pino, Kowalczyk and Wei in [21] that shows the dimension eight is the critical dimension. For 
dimensions 4 < n < 8 there are various partial results under certain extra (natural) assumptions on solutions 
by Ghoussoub and Gui in [37], Savin in [47] and references therein. The remarkable point is that in lower 
dimensions, that is when n < 3, this conjecture holds for the scalar equation 

(1.12) —A u = H(u) in R", 

with a general nonlinearity H G C ll (R). In this regard, we refer interested readers to [36] by Ghoussoub and 
Gui for n = 2 and to [2] by Alberti, Ambrosio and Cabre for n = 3 and also to [28] by Farina, Sciunzi and 
Valdinoci for a geometric approach in two dimensions. When the Laplacian operator in (1.12) is replaced 
with the fractional Laplacian operator (—A) s , De Giorgi type results are known when either n = 2 and 
0<s<lorn = 3 and 1/2 < s < 1, see [8-10, 49]. However, in three dimensions and when 0 < s < 1/2 the 
problem seems to be more challenging and remains open. 

More recently, Ghoussoub and the author in [31] provided De Giorgi type results up to three dimensions 
for the following elliptic gradient systems when u : R" —»• R m 

(1.13) -A u = VH(u) in R”, 
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for a general nonlinearity H £ C 2 (R m ). In this regard, authors introduced the notions of orientable systems 
and //-monotone solutions to adjust and to apply the mathematical techniques and ideas given for the scalar 
equation (1.12) to (1.13). It seems that these concepts are essential, in this context, to explore system of 
equations. Note also that Sire and the author in [32] studied system (1.13), when the Laplacian operator 
is replaced with the fractional Laplacian operator (—A) s , and provided De Giorgi type results for certain 
parameters s = (si, • • • , s n ) £ R” in lower dimensions. 

Definition 1 . 2 . System (1-9) is called orientable, if there exist nonzero functions 9k £ C 1 (R™ +1 ), k = 
1, • • • ,m, which do not change sign, such that for all i, j with 1 < i < j < m, 

(1.14) djHi{u)0i(x)0j(x) > 0 for all x £ 1”. 

Similarly, if the condition (1-14) holds for the extension v, then we say (1-9) is orientable. 

Note that the orientability imposes a combinatorial assumption on the sign of the nonlinearity H = (Hi)” 1 1 
and therefore on the system. 

Definition 1.3. We say that a solution u = (ui)7S 1 of (1.1) is H-monotone if the following hold, 

(1) For every i £ {1, • • • , m}, Ui is strictly monotone in the x n -variable (i.e., d n Ui ^ 0). 

(2) For i < j, we have 

(1.15) djHi(u)d n Ui(x)d n Uj(x) > 0 for all x £ R". 

Similarly, if above conditions hold for the extension function v then we say that v is H-monotone. 

It is straightforward to observe that /Z-nronotonicity implies stability. One can show this by differentiating 
(1.1) with respect to the variable x n and setting <fi := d n Ui. As the last definition of this section, we provide 
the notion of symmetric systems that plays a key role in our main results. 

Definition 1.4. We call system (1.1) symmetric if the matrix of gradient of all components of H that is 

(1.16) H := (djHiiu))^ 

is symmetric. Similarly, if H := (djHi(v))™j =1 is symmetric for the extension v, then we say (1.9) is 
symmetric. 


In this article, we prove De Giorgi type results for stable and //-monotone solutions of symmetric system 
(1.9) with a general nonlinearity H when n = 2 and 0<Sj<lorn = 3 and 1/2 < s* < 1, see Theorem 
4.2. We also provide Liouville theorems for stable solutions of (1.9) under some extra assumptions on the 
nonlinearity H when n < 2(1 + s*), see Theorem 4.4. In addition, we establish the following geometric 
Poincare inequality for bounded stable solutions of (1.9). Let r\ = (ilk)™=i £ C/(R” +1 ), then 

(1-17) it-T [ y 1 - 2 ^ (|Vuj| 2 ^lf + |VT«|V x t>j|| 2 ) rftdx 

i= i /{|v a ,w|#o}nR" +1 

(1.18) -fV' / ( J djHi(v)diHj(v)\\7 x Vi\\V x Vj\r]i7y - </,//,( r)Y , r, • V x Vj-qf ) dx 

Jd K +1 V J 


s HA 


\V x Vif\S/rn\ 2 dx, 


where Vr, stands for the tangential gradient along a given level set of Vi and .4 2 for the sum of the squares 
of the principal curvatures of such a level set, see Theorem 3.2. Note that for the case of scalar equations, 
m = 1, a similar inequality was established by Sternberg and Zumbrun [50] to study phase transitions and 
area-minimizing surfaces. For this case, the boundary term (1.18) disappears. The idea of applying the 
geometric Poincare inequality to prove De Giorgi type results was initiated by Farina, Sciunzi, Valdinoci in 
[28] and references therein. Ghoussoub and the author in [31] established a counterpart of this inequality 
for local systems of the form (1.13). Note that for symmetric systems, m > 2, the boundary term (1.18) 
becomes 


( 1 . 20 ) 


y / djHi(v) (\V x Vi\\W x Vj\rnrjj -V x Vi ■VxVj'nf) dx, 
771 Jo R " +1 
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where the integrand has a fixed sign for an appropriate test function ry = (rji)^_ l . In the light of this 
inequality, we prove De Giorgi type results in two dimensions and we show that vectors \7 x Vi(x, 0) and 
V x Vj(x, 0) for i ^ j are parallel and the angle between two vectors is arccos ( ) • 

Consider the scalar equation (1.12) when H > 0 for H'{u) = —H(u). Modica in [44] proved that the 
following pointwise estimate holds for bounded solutions 

(1.21) |Vw| 2 < 2 H(u) in R”. 


This inequality has been used in the literature to study entire solutions of semilinear elliptic equations and, 
in particular, to establish De Giorgi type results, see [1, 2, 4, 21, 36, 37, 39, 47]. Unfortunately, a counterpart 
of this inequality does not hold for the system of equations of the form (1.13) with a general nonlinearity. 
However, assuming that n = 1 and m> 1 and multiplying the i th equation of (1.13) with v! i and integrating, 
it is straightforward to observe that the following Hamiltonian identity holds 

1 m 

(1.22) 2^l M ’( a ')l 2 +H(u(x))=C for irt, 


when C is a constant. Gui in [39] examined system (1.13) when n > 1 and m > 1 and proved the following 
elegant Hamiltonian identity in higher dimensions 


(1.23) 


m ^ 

9 (|V*'Ui( a 


\d n Ui(x)\ 2 ) - H(u(x)) 


dx' = C for x n € R, 


where x = ( x',x n ) € R ra and C is a constant. In this paper, as Theorem 3.1, we provide a counterpart of 

(1.23) for the fractional system (1.9) in one dimension that is 


(1.24) 



( dyVi ) 2 dy — H(v(x, 0)) = C for xel, 


when C is a constant and diH{v) = Hi(v) for every i. Proving a similar identity in higher dimensions, n > 2, 
remains an open problem. Our methods and ideas are strongly motived by the ones provided by Gui in [39], 
Cabre and Sola-Morales in [10], Cabre and Sire in [8, 9] and Sire and the author in [32]. 

Lastly, we consider a two-component system of the form (1.1) with the following particular nonlinearity 

(1.25) Hi(m, u 2 ) = H\u\ + j 3 u\u\ and H 2 (ui, u 2 ) = H2u\ + Pu\u2- 

where u : R" —> R 2 and yi,H2,P are constants. Both local and nonlocal systems of this type, known as 
the nonlinear Schrodinger system, have been studied extensively in the literature. For more information, we 
refer interested readers to [3, 15, 29, 43, 46, 51-54] and references therein. The nonlinear Schrodinger system 
is a natural counterpart of the following nonlinear Schrodinger equation, 

(1.26) Au — u + u 3 = 0 in R". 

Even though equations (1-26) and (1.11) look alike, their solutions behave very differently. In this article, 
we provide various Liouville theorems and nronotonicity formulas for solutions of the nonlinear fractional 
Schrodinger system, under some assumptions on parameters Si, S 2 , Hi, 1 ^ 2 , P- 

The organization of the paper is as follows. In the next section, we provide some standard regularity 
results and estimates regarding fractional Laplacian operator. In Section 3, we prove various technical tools 
needed to establish our main results. We provide a Hamiltonian identity, a geometric Poincare inequality 
and various gradient and energy estimates for system (1.9) with a general nonlinearity H = (iT,)™_ 1 . We 
also provide a monotonicity formula to conclude the optimality of gradient and energy estimates. In Section 
4, we establish De Giorgi type results and Liouville theorems for symmetric systems in lower dimensions. In 
addition, applying the Hamiltonian identity we analyze directional derivatives of the extension function v 
satisfying (1.9). Lastly, in Section 5, we consider a two-component fractional Schrodinger system, that is a 
particular case of (1.1), and we provide various Liouville theorems and monotonicity formulae. 
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2. Standard Elliptic Estimates 


In this section, we provide some standard estimates regarding fractional Laplacian operator. We assume 
that the nonlinearity H = {Hi)" i 1 and each Hi belongs to C' 1 ’ 7 (R m ) with 7 > max(0,1 — 2s*) and u = ( Ui) r [L 1 
is a sequence of bounded functions. We omit the proofs of lemmata in this section and we refer interested 
readers to [6-9, 32, 49] and references therein. We start with the following regularity result for solutions of 
(1.1). 

Lemma 2.1. Suppose that u = ( uf)i is a bounded solution of (1.1). Then, each Ui is C ,2, ^(K") for some 
0 < /? < 1 depending only on Si and 7. 

For the sake of convenience, we use the following notation used frequently in the literature. 

Notation 2.1. Set B R = {x = (x, y ) G R” +1 , |x| < R}, d + B R = 8B R D {y > 0}, V° R = dB R n {y = 0} and 
C H = B n x (0. R) for Ii > 0. 

We now provide a regularity result for solutions of an equation in the half-space. This can be applied to 
extensions functions v satisfying (1.9). 

Lemma 2.2. Let fi G C^IT^r) for some a G (0,1), R > 0 and Vi G L°°(B 2 r) D H 1 {B^r, y 1-28 *) for some 
Si G (0,1) be a weak solution of 

/ di viy^Vvi) = 0 in B+ R cWf + \ 

{ -lirn^o y l ~ 2si ^ = fi on r° fl . 

Then, there exists 0 G (0,1) depending only on n, Si, and a, such that Vi G C 0 , ^(B R ) and z/ 1 2si <9^^ G 
C°'P{B r ). Furthermore, there exist constants Cr and Dr such that 

( 2 - 2 ) - Ci ^ and W yl ~ 2Sid y Vi \\co,f>(B+) - Di ( R )’ 

where Ci{R) only depends on n,ai,R,\\vi\\L°°{B2 R ) and ||/i||z^°°(rg_ R ) an d Di(R ) oriZz/ depends on n, at, R, 
\\vi\\ L oo {B+ r ) and\\fi\\c°(ro R) . 

The next lemma provides gradient estimates for bounded solutions of (1.9) and it is a consequence of 
Lemma 2.1 and lemma 2.2. For more information, we refer interested readers to [7-9, 48] and references 
therein. We apply these estimates frequently in the proofs of our main results. 

Lemma 2.3. Let v = {vi) r fL 1 be a bounded solution of (1.9). Then, each Vi satisfies 


(2.3) 

<1 

8 

c* 

IT 

< 

C for x G I” and y > 0, 

(2.4) 

|Vuj(a;, y)\ 

< 

. , for (x,y) G R" +1 , 

1 + y 

(2.5) 

\y l ~ 2s id y Vi(x, y)\ 

< 

C for (x, y) G R" +1 , 


where the positive constant C is independent from x and y. 


3. Hamiltonian Identity and Analytic and Geometric Estimates 


We start this section with proving a Hamiltonian identity for solutions of (1.9) in one dimension. We 
then apply this identity, see Theorem 4.4, to study bounded stable solutions of symmetric (1.9) with a 
general nonlinearity H = (7/, t )™ :1 . Hamiltonian identities are amongst the most important tools to study 
qualitative behaviour of entire solutions of differential equations and in some cases they lead to properties 
such as monotonicity formulae. 


Theorem 3.1. Suppose that v = {vi) r jL 1 is a solution of (T9) in one dimension and 0 < s.j < 1 for 
1 < i < w. Then the following Hamiltonian identity holds 


(3.1) 



( dyVi ) 2 dy — H(v(x, 0)) = C for iGl, 


where C is a constant that is independent from x and diH(v) = Hi(v ) for every i. 
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Proof: Let v = (t>j)™ 1 be a solution of the extension problem (1.9). Set 


m -i ,00 

'j( x ) : = J o y 12si K d ^ v i ) 2 - ( d v v i) 2 ] d v- 


Differentiating w in terms of x, we get 


m 1 f°° 

d x w(x) := V] — / y l ~ 2si [- d x Vid xx Vi - dyVid xy Vi] dy. 

: , JO 


From (1.9) for each i = 1, • • • , ?n, we have 

(3.4) y l ~ 2si d xx Vi + (y 1_2si <9 y u,:) = 0. 

Combining this and (3.3), we end up with 


171 1 ,00 

d x w(x) := -T- / (y 1-284 ^*) - dy. 

;_I d Si J 0 


Integration by parts for the first term in the above yields 


,oo ,oo 

/ d x Vidy (y 1 - 28 * dyVi) dy = / y 1-2 ^d xy Vid y Vidy - lim y l ~ 2si d y Vid x Vi. 

Jo Jo y^° 


Substituting this in (3.5) we get 


d x w(x) = - — lim y 1 2si d x Vid y Vi. 

^ d Si ?/-*■ o w 


From this and the boundary term in (1.9), we have 


Therefore, 


m 

d x w(x) = ^2diH(v(x, 0))d x Vi = d x [h{v{x, 0))) . 

i-1 

d x w(x) — H(v(x, 0)) =0. 


This implies that w(x) — H(v(x, 0)) is constant in terms of x. 

□ 

We now prove an inequality, known as the stability inequality, for stable solutions of (1-9). This inequality 
plays an important role in this paper. 

Lemma 3.1. Let v = {vi) 1 fL 1 be a stable solution of system (1.9). Then, the following stability inequality 
holds 

m ,, _ m ,, 

(3.10) "Y / Jd Si d Sj djHi[v)d Vi Hj(v)QQdx <Y I y ai \\/Q\ 2 dx, 

i,j =l Jd K +1 i=l J K +1 

for all Q G C 1 c (R n + + 1 )- 

Proof: Since v = (ty)”1 1 is a stable solution of system (1.9), there exists a sequence (f> = that satisfies 

c 2 

(1.10) . Multiply the i th equation of (1.10) with and do integration by parts to get 

(3.11) - / y° 4 V<^ ■v( < J-)dx+ [ y ai V<j>i ■v(^-)dx = 0. 

JrI +i KfaJ Jd r" +1 

From this and the boundary term of equation (1.10), we get 

(3.12) -2 / y° 4 V^-VCi$-dx+ / y ai |V^| 2 -Jda;+ / d Si Y djHi{v)<t>Adx = 0. 

J R^ +1 <Pi JRI+ 1 <Pi JdRl +1 “ <Pi 
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Applying the Young’s inequality, for each index i, we obtain 


( 3 . 13 ) 


[ d Si y djHiivfetfdx = f y ai + 2 V& • VCiT-) 

Jd r ” +1 (pi J R "+! V <pf <pij 

< f y ai \V(i\ 2 dx. 

J R " +1 


We now provide a lower bound for the integrand in the left-hand side of (3.13) when taking the sum on both 
sides of (3.13) for i = 1, • • • , m, 

(3.14) yd Si djHi{v)^ - ydsAH^vy+yds^H^v)^ 

hJ * 

= y dsAHiiv) C 2 + 5] d 8i djHi(v)^-Ci + y dsAHA)^ 

i i<j ( P‘ l i>j V 1 

= y dsA H iiv)(i + y d Si dj H i(v)^-C 2 + y d Sj diHj(v)^C 2 

*Pi , 

l 1 <J 1 <J J 

= y dsA H i( v Ki + y^itpjY 1 {dsidj H i(y)<t> 2 j(,i + d aj diHj{v)^Cj) 

i i<j 

> yd s AHA) c 2 + 2E^ ds. dj II, (v)diHj (v)QQ 
i i<j 

= y \Jd Si d Sj djHAAH]Y)CiC,o- 

i,3 

This completes the proof. 

□ 

We now apply the stability inequality (3.10) to establish a geometric Poincare inequality. Note that this 
inequality for the case of local scalar equations was first proved in [50] and it was used in [28] and references 
therein to prove De Giorgi type results. For the fractional Laplacian case this inequality was established 
in [49]. We also refer interested readers to [31] and [22, 32] for the case of local and nonlocal systems, 
respectively. 

Theorem 3.2. Assume that m, n > 1 and v = (iq)i is a stable solution of (1.9). Then, for any y = 
(rf = i £ C^(R" +1 ), the following inequality holds, 


i=1 ddv^il^oinR" 4 


y 1 2si (|Vu,| 2 . 4 2 + IVtJV^H 2 ) y 2 dx 


Jd K +1 


idjHi(v)diHj(v)\V x Vi\\V x Vj\r]ir]j - djHi(v)X7 x Vi ■ X7 x VjT / 2 ) dx 


m i r 

< V— / y 1 - 2 si |V x t) i | 2 |Vr?j| 2 dx, 

“ d Si Jr ” +1 


i=l 


where V t, stands for the tangential gradient along a given level set of Vi and A 2 for the sum of the squares 
of the principal curvatures of such a level set. 

Proof: Let v = (vi) r fL l be a stable solution of ( 1 . 9 ). From Lemma 3 . 1 , the stability inequality ( 3 . 10 ) holds. 
Test the stability inequality with Q := \V x Vi\r]i where each y 1 £ C,i(M” + ), to get 


' == E 


' d Si d Sj djHi(v)d Vi Hj (v) \ X7 x Vi \ \ X7 x Vj\yiyjdx 


i,j =1 “ 


m !> 

< y y ai \X7 {\X7 x Vi\rn) \ 2 dx =: J. 
, J R" +1 
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Simplifying the left-hand side of the above inequality, we get 


m « 

I = E / d Si d Vi Hi(v)\W x Vi\ 2 7^dx 

Yl J dK +1 


JdK + +1 


' d Si d Sj djHi ( v)d Vi Hj (u)| V x Uj| | V x Vj\r)iT)jdx. 


Similarly, for J we have 


m r r 1 

(3.18) J = Y y ai \^xV i \ 2 \V Vi \ 2 + \V\V x v i \\ 2 r 1 f + -V\V x v i \ 2 -Vr ] f dx. 

j_ 1 JR ” +1 L ^ 

We now differentiate the i th equation of (1.9) with respect to Xk and multiply with d Xk Virff for i = 1, • • • ,m 
and k = 1, • • • , n, to get 

(3.19) div(yfV d Xk Vi)d Xk Virjf = 0. 

Integrating this by parts, we obtain 


[ y at \^dx k v i \ 2 ri 2 i dx + \ [ y ai V\d Xk Vi\ 2 • Vrftdx 

y R n + i l y R n + l 

/ lim y ai (- dyd Xk Vi)d Xk v^dx. 

Jd r " +1 y^° 


Differentiating the boundary term of (1.9) with respect to Xk yields 


(3.21) 


lim y ai {-d y d Xk v i )d Xk v l ri 2 i = d Si Y d J H i (v)d Xk v j d Xk v l r) 2 . 

3 =1 


From this and (3.20), we have 


m m ^ p 

Y / y ai \Vd Xk Vi\ 2 rj?dx + Yo y ai ^\d Xk Vi\ 2 ■ Vrjfdx 

J K +1 t ! 2 M +1 

m p m p 

= Y d Si d Vi Hi(v)\d Xk Vi\ 2 r]fdx + Y d Si / djHi(v)d Xk Vjd Xk Vir] 2 dx. 

Jd R " +1 Jo R ” +1 


<=1 ^R" +1 


JdR + +1 


Taking sum on the index k = 1, • • • , n, we get 


m p m p n 

(3.23) Y ds * / d Vi Hi(v)\V x Vi\ 2 rijdx = Y V a> E \Vd Xk v i \ 2 rj 2 dx 

i= 1 Jm l +1 i= i J K +1 k=i 

1 m r 

+oE / y^viv^i 2 • Vrftdx 

2 M +1 


-Y d ° 


i#* JaR + +1 


djHi(v)X7 x Vi ■ V x Vjr] 2 dx. 


We now substitute this equality in I < J when I is given in (3.17) and J is given in (3.18). On can see that 
the term \ YILi / R n+i y ai VlV^I 2 • Vrjfdx crosses out and we end up with 


(3.24) Y y ai E - |V|V^|| 2 y 2 dx 

i= 1 - / »? +1 n{|Vt>,| ? 40} \fc =1 J 


(3.25) +E / \/d Si d Sj d v H i (v)d Vi H j (v)\W x v i \\W x v j \r] i r] j -d Si d j H i (v)W x Vi-'V x v j ri 2 dx 

U Jd R + 

m p 

(3.26) < Y / 
i =i d K +1 


y ai \V x Vi\ 2 \\/m\ 2 dx. 



According to formula (2.1) given in [50], the following geometric identity between the tangential gradients 
and curvatures holds. For any w £ C 2 (fl) 


(3.27) ^ \Vd k w\ 2 - |V|Vw|| 2 

k= 1 


n 2 ) + |Vt|Viu|| 2 for x £ {|Vru| > 0 D fl}, 
0 for x £ {|Vu>| = 0 fl fi}, 


where Ki are the principal curvatures of the level set of w at x and Vt denotes the orthogonal projection of 
the gradient along this level set. Applying the above identity (3.27) to (3.24) completes the proof. 

□ 

For the rest of this section, we provide energy and gradient estimates for solutions of (1.9). Then, in next 
sections we apply these estimates to establish De Giorgi type results and Liouville theorems. Consider the 
energy functional 

yi _2si |Vr’i| 2 dx — f H(v)dx, 

Jb r 


(3.28) 


«»> '=Y.Ti~ 


>C 


when diH(v ) = Hi{v ) for i = 1, • • • m. We finish this section by proving an energy estimate for //-monotone 
solutions of (1.9). 


Theorem 3.3. Let v = (vif^L 1 be a bounded H-monotone solution of (1.9) such that 

(3.29) lim Vi(x',x n ,y) = Li for i' £ l"' 1 and y £ R + , 

x n —>oo 

when L = (Li) r fL 1 and each Li is a constant in R and H(L) = 0. Then, the following energy estimates hold. 

(i) 1/1/2 < s* < l, then E R (v) < CR n ~ 1 , 

(ii) If s* = 1/2, then E R (v) < CR ™ _1 log R, 

(iii) J/0 < s* < 1/2, then E R (v) < CR n ~ 2s *, 

where the positive constant C is independent from R > 1. 

Proof: Set the shift function uf (:r, y) := Vi(x', x n +t, y) for (x 1 , x n ,y) £ R" +1 and t £ R. It is straightforward 
to see that v 4 = is a solution of (1.9) and it satisfies pointwise estimates provided in Lemma 2.3. In 

addition, for every parameter t and all indices i, one can see that |i>|| £ L°°(R/ +1 ). Therefore, v * = (i>f)’ji 1 
is a sequence of bounded functions and 

(3.30) lim \v*(x,y) — LA + \Vvl(x,y)\ = 0 for all (x, y) £ R+ +1 . 

£->-oo 

This implies that 

(3.31) lim E R (v*) = 0. 

£—>• oo 

We now differentiate the energy functional in terms of parameter t to get 

^ 1 pR P p TYl 

(3.32) dtEn^) = 53—/ / y^VvCVidtvDdx- 

® s i JO JBr j Br 

m -| pR p 

= ~Y.fr / div (y^ 2s ‘Vu‘) d t v\dx 

i—l Q'Si JO J Br 
rrL i pR p 

+ E T / / y'~ 2Si ( Vu J • v ) dtvldH^dy 

® s i JO JdB R 

m „ 

-53 / H i (v t )d t vjdx. 
i=1 Jb r 

From the fact that div (y 1_2s *Vu|) = 0, we can simplify the above as 

-j pR p etl /* 

(3.33) d t E R (v t ) = 'S2— / y 1 ~ 2si 'Vvi-udtV t i d'H n ~ 1 dy + 'S2— y^^dyvldtvjdx. 

l = l ® JS i JO JOBr dsi JB R x{y—R} 
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Consider disjoint sets of indices / and J such that / U J = {1, • • • , to} and d t v * > 0 > <9iU A for y G I and 
A € J. Applying this, we can expand dtE^v*) as 


(3.34) &£«(«*) = / / E + E dydn 

JdBn Jo \ a s u ,, «s A / 


E -^-y 1 - 28 ^*^* + E j-y 1_28A ^i/UA^U ate. 


JS B x{y=it} y Me/ u V ASJ SA / 

Lemma 2.3 implies that there exists a constant M such that \d v v\\ < for (x, y) € R" +1 and |9 y u|| < 
for x G 1" and y > 1. From these estimates and (3.34), we get 

(3.35) dtEniv*) > M f V -^y l ~ 2s » f-E— ) dydU^ 1 

JdB R J o V n-yy M 


mg/ 

!>R 1 


EA 


/obr Jo Ag , u s x 




d t v\dyd'H 


Ib r x{ v =r} m6J “s 


Er» M ' 


Er»“‘ OWl*- 


JB H x{y=B} A£J \y/ 

Note that E r (v) = E R (v T ) — / 0 T dtE^v^dt for every T > 0. Combining this and (3.35), we conclude 


Er(v) < E R (v T ) + M 


f [ R y 1 ( y 1 ^ 
ldB R Jo fJ , eI ds ^\ 1 + y 


dtV^dtdydH 7 


1 / \ [ T 


ldB R Jo X( -J dsx V 1 + y 


d t v x dtdydH T 


E J-y 1 ' 28 ^ 1 I dtv^dtdx 


IB R x{y=R} ^ ds 


Simplifying (3.36) yields 


/ E— y 

/fl«x{9=fi) X& J ds x 


P pR ]__ 2 s 

(3.37) £*(«) < E r (v t ) + M / ^-—^{vl-v^dydW 

JdB R J o “ ! + y 

/• „ ,.i- 2 s a 

+M / E - v x )dyd'H n ~ 1 

Job r Jo r 1 + y 


1 „.l-2s A -l 


d t v x dtdx. 


+M l E y 2sm _ v ^) dx + M [ E y 2s * ( wa - 

JB K x{j/=i?} JSrx{i/=K} Ag j 

Note that for p € I and A € J, we have v x > v M and v\> v x . Therefore, 


(3.38) £i?M<£^ T ) + M 


/* /*-R m _ —2si /* 

/ / E Y—dydH n ~ 1 + M Y,y~ 2Sidx - 

ldB R Jo 1 + y Jb r x{ V =r}^ r[ 
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Let T —>■ oo, then E R {v T ) approaches zero. Doing integration by parts we obtain 

m 

(3.39) E r (v) <MJ2 [R n ~ 2si X{o< Si <i/ 2 } + R n ~ 1 X{ i/ 2 < Si <i } + R n ~ x logi? X{s ^i/ 2 } + R n ~ 2si ] ■ 

i=\ 

This completes the proof. 

□ 

Note that ideas and techniques applied in the above proof are strongly motivated by the ones provided 
by Ambrosia and Cabre in [2] and by Cabre and Cinti in [6, 7] for local and nonlocal scalar equations, 
respectively. We would like to refer interested readers to [31] by Ghoussoub and the author and to [32] by 
Sire and the author for local and nonlocal system of equations, respectively. In the next two lemmata, we 
prove gradient estimates for bounded solutions of (1.9). 

Lemma 3.2. Suppose that v = (fi)i^i a bounded solution of (1.9). Then, the following gradient estimate 
holds for any ?' = !,••• , to 


(3.40) 


[ y ai \X7vi\ 2 dx<CR n 
JC R 


when C is a positive constant independent from R. 

Proof: Let v be a bounded solution of (1.9). Multiplying (1.9) by Viip R where i/jr is a test function and 
integrating over K" + , we get 


(3.41) 


y a *\Vvi\ 2 ipRdx + 2 y ai SJi>R ■ '\/v i ip R v i dx = d s 


Hi(v)viip R dx. 


Using the Cauchy-Schwarz inequality, we conclude 


(3.42) 


y ai if R Vil)R ■ VviVidx < - 

>i+i 2 


y ai \Vui\ 2 if 2 R dx + 2 / y ai \S7ip R \ 2 vfdx. 


From this, we obtain 

(3.43) f y ai \Vvi\ 2 ip%dx < 4 


y ai \S7ip R \ I v(dx + 2d. 


/9R7 


Hi(v)viip R dx. 


From the boundedness of v and the fact that Hi £ C 1 ’ 7 , we obtain the following bound by choosing an 
appropriate test function, 


(3.44) 


>c R 


y ai \^7vi\ 2 i/j R dx < CR 


-yn—2si 


+ CR n < CR n , 


where the constant C is independent from R but may depend onw= (vi) r (L 1 . This completes the proof. 

□ 

We now assume extra conditions on the sign of the nonlinearity H. This enables us to prove stronger 
gradient estimates on solutions of (1.9). 

Lemma 3.3. Suppose that v = (r;,)^. 1 is a bounded solution of (1.9). 

(i) If Hi(v) > 0 for any 1 < i < m, then J Cr y 1 ~ 2si \Vvi\ 2 dx < CR n ~ 2si . 

(ii) IfY!iLi v i H i( v ) < then YhLi I Cr y l -' 2si \Vvi\ 2 dx < CR n ~ 2s -. 

Here, the positive constant C is independent from R. 

Proof: We start the proof with Part (i). Suppose that for a particular index 1 < i < m we have Hi(v) > 0. 
We now multiply (1.9) with (vt — ||r , i||oo)V’fl when ijj R is a positive test function. Applying integration by 
parts and using the fact that Hi has a fixed sign, we obtain 


(3.45) / y ai 'V[(vi-\\vi\\ 00 )ip R }-'Vvidx = / \mny ai (-d y Vi){vi-\\vi\\ 00 )'ipRdx 

JR" +1 J0R” +1 y^° 

r 

| \ ^i 11 oo )Hi(v)ip R dx < 0. 


= d s 


n 



This implies that 


(3.46) [ y ai \Vvi\ 2 ip R dx < 

J R” +1 


- / y ai {Vi - ■ WVidx 

J R” +1 

- I 2/ ai V(ui - Halloo) 2 • 

2 J R” +1 

f {Vi -\\viWoo) 2 dw(y a ‘Vi/j R )dx 
4 JR" +1 

+ 7: / (- lim y ai dyip R )(vi - ||^||oo) 2 <ia: 

4 Jd R” +1 J /_>0 


—: I + J. 

From boundedness of Vi and applying an appropriate test function ip R , we obtain 
(3.47) I = -- 

p2R / n.Cii n.Cli — 1 


4/ K - Halloo) 2 {y ai A x ip R + aiU ai 1 d y ip R )dx 
Z jR n+1 


'.. c 


' C 2 R 


y ai y ai ~ l 
R 2 R 


dx < CR n 


IR 


y 1 , y 

R 2 


R 


dy < CR 


n—2si 


when the positive constant C is independent from R. In above, we have used the fact that the test function 
ip R can be chosen such that \d y tp R \ < CR~ X and \A x ip R \ < CR ~ 2 in C R . Using these properties, we conclude 
\y ai dyil ) R | < CR~ 2si in C R . Similarly, applying the test function \p R with the above properties, we obtain 
the following upper bound for J given in (3.47) 

(3.48) J < CR n ~ 2si . 

This completes the proof of Part (i). We now provide a proof for Part (ii). Multiplying the i th equation of 
(1.9) with ViiJj R , when ip R is the same test function applied in Part (i), and integrating we get 


(3.49) 


m 1 « m ^ 

V'—/ y ai V[viip R }-Vuidx = V' / lim y at (-d y v i )v i 'ip R dx 
d Si J R»+! Jd R»+! y^o 

P m 

= / T ViHi(v)ijj R dx. 

JdRl +1 i=1 

Note that the latter is nonpositive since BhLi v iHi( v ) < 0. The rest of the proof is similar to Part (i) and 
we omit it here. 

□ 

We end this section with a nronotonicity formula for bounded solutions of (1.9). Consider the following 
function for R > 1, 


(3.50) 


I{R) := R 


2 s* —n 


1 


yy ai \X7vi\ 2 dxdy- 


3 + 

^R i= 1 


1 B R x{y= 0} 


H{v)dx 


when diH{v) = H-i(v) for every i. We show that I(R) is a nondecreasing function of R, when H < 0, that 
is equivalent to 

(3.51) I(R) > 7(1) for R > 1. 

Note that from definitions of I(R) in (3.50) and E R (v) in (3.28), we have 

(3.52) E r (v) > i? n - 2s *J(i?). 

From (3.52) and (3.51), we conclude 

m 1 r r 

(3.53) y ■ 


i=1 2d Si Jc R 


y l ~ 2Si \Vvi\ 2 dx - / H(v)dx > CR 


n—2s* 


’ Br 


and the positive constant C is independent from R. This clarifies the sharpness of gradient and energy 
estimates provided in this section, in particular when 0<s* = s*<l/2. 
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Theorem 3.4. Let v = (fi)i2=i be a bounded solution of (1.9). Then, I(R ) is a nondecreasing function of 
R> 1 when diH(y) = Hi(v) for every 1 < i < m and H < 0. 


Proof: Differentiating I{R) with respect to R, yields 
n - 2s* 

iBi 


(3.54) I'(R) = - n 0 2 ' S Ji 2s *- n ~ 1 f V'' y ai \Wvi\ 2 dxdy + (n — 2s*)R 2s *~ n ~ 1 f 

2 Jb± , n JE 


SRx{y=0} 


H(v)da 


1 n m ,. 

+ R 2s*-n S"y ai \S7ui\ 2 dxdy-R 2s ~ n / //(t-)(ix. 

2 Jd+B+^ JdB R x{y= 0} 

Applying some standard arguments in regards to the Pohozaev identity implies that 


(3.55) 


, t 7Tt 7Tb a /» a 

5/ Vy ai | y ^| 2 = Vs., / y°‘|Vt>i| 2 +n / H(v)dx-R H(v)d'H n ~ 1 

Z J B+ , , , =1 dB+ JBr 


_ -R / E y 0i {d vVi ) 2 dH n + f / E^* I 

" 2 ^+s+ i=1 


7 9+^ i=l 


Combining (3.54) and (3.55), we conclude 

« m 

(3.56) /'(R)lT +1 - 2s * =R E y ai (dvV i ) 2 d'U n 

Jd+B+ ' 


« 771 r. 

+ y^is* - Si )y ai \Vvi \ 2 - 2s* H(v)dx. 

dB~f d B R x{y— 0} 


Since H < 0 and s* > s,; for every 1 < i < m, we have /'(l?) > 0. 


□ 


4. Symmetry Results and Liouville Theorems; General Nonlinearity 


In this section, we provide De Giorgi type results for stable and fJ-monotone solutions of symmetric 
system (1.9) with a general nonlinearity in lower dimensions. Just like in the proof of the classical De Giorgi’s 
conjecture in dimensions n = 2,3 providing a Liouville theorem for the quotient of partial derivatives is a 
milestone. To be mathematically more precise, consider the case of scalar equation (1.12), it was observed 
by Berestycki, Caffarelli and Nirenberg in [4], by Ghoussoub and Gui in [36] and by Ambrosio and Cabre in 
[2] that if <f) £ Lj oc (R") such that fa 2 > 0 a.e. and a £ Hf oc ( R") satisfy 

(4.1) —crdiv(0 2 cr) < 0 in 1", 


under the decay-growth assumption 


(4.2) / </>V < CR 2 for R > 1, 

J B2r\Br 

then a must be a constant. Eventually, a is set to be a = EfT f° r an arbitrary direction r\ £ R” and 
<f> = d n u. Cabre and Sire in [8, 9] proved a similar Liouville theorem for nonlocal scalar equations that is 
(1.1) for m = 1. Most recently, Ghousssoub and the author in [31] and later Sire and the author in [32] 
provided counterparts of this Liouville theorem for the local and nonlocal gradient system (1.13), respectively. 
Consider the following set of functions T that is somewhat standard in this context 

f°° dr 

: R + —> R + , F is nondecreasing and / -— = c 

J 2 rF(r) 

For example, F(r) = lnr and F(r) = constant > 0 belong to this class and F(r) = r does not belong to J. 
As far as we know, this set of functions was introduced by Karp in [40, 41] and used in [6-9, 30, 32, 45]. 


F:= <F 


Theorem 4.1. Assume that fa £ L[^ C (M" +1 ) is a positive function and a £ ,y ai ) satisfies 


(4.3) 


lim sup 

Ft—yoo 


i 

R 2 F(R ) 


m 

< °°> 
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(4.4) 


for F £ T and i = 1, • • • , to. If the sequence (ai)'fL 1 is a solution of 

—<7; div(y ai (f>?V(Ji) < 0 in R" +1 , 

-lim^o y ai <f\(Jidy<Ji < Z)" =1 - a i) a i in 3M” +1 , 

when 0 < hij £ Lj oc (M. n ), hij = hj t i and f £ Lj oc (M.) is an odd function such that f(t) > 0 for t £ R + . 
Then, each function <x; is constant for all i = 1, • • • , to. 

In this article, we apply the above theorem frequently to prove our main results for solutions of (1.9). 
Suppose that v = (r^ELi is a //-monotone solution of (1.9). Let fa := d n Vi and ifi := V'c, • rj for any fixed 
r] = (f/,0) £ R" _1 x {0}. Then, sequences ( wf) r fL 1 = (fa)il i and {wi) r fL 1 = ('ifi) r fL 1 satisfy the following 
linearized equation 

div(y ai Vuy) = 0 in R" +1 , 

-lim^o y ai d y Wi = d Si J2™=idjHi(v)wj in <9R” +1 . 

Since v is a //-monotone solution, fa does not change sign for each i. Then, the sequence of functions 
cr = (o-j)i(Li f° r a i : = ^7 satisfies the following 

(4.6) di v(y ai 4qVai) = div(y 0i [faVai - a,Vfa]) = fa di v{y ai Vfa) - fa di v(y a 'Vfa) = 0. 

On the other hand, in <9R(|_ +1 we have 


(4.5) 


(4.7) 

which implies that 


I m 

lim y ai dyOi = - lim y^dyfa^f 1 + lim y ai d v fa~2 = d Si Y' djH,{v) 
y-> 0 y-> 0 y-> 0 07 f 

1 3 —1 


fa 


, i’i 


(4.8) - lim y ai fa\<jidy(Ji = d Si 22 djHi(v)[fafa - (fjfafa = d Si 22 </,//, (r)0;0,i' T , - a l )a i . 

3= 1 3 -1 

Note that if we set faj = djHi(v)faf)j and / to be the identity function, then the above equations (4.6) and 
(4.8) satisfy (4.4). Note that for symmetric systems we have hij = hj t i and for //-monotone solutions we 
have hi,j > 0. The following computation for the right-hand side of (4.8) is our main observation to define 
symmetric systems and //-monotone solutions and to establish Theorem 4.1; 


(4.9) 22 fa <A? 9jHi(v)ai(ai - crj) = 22 fafa<djHifafatfai - aj) + 22 fafadjH t (v)ai(ai - aj) 

i,j=l i<j i>j 

= 22 fafpjdjH^faaifa - + 22 OiOjiflfiOa, ( a 3 - af) 

i<j i<3 

= '22 fafadjHi(v)(ai - aj) 2 > 0 . 

i<j 

We are now ready to state the following De Giorgi type result. 

Theorem 4.2. Suppose that v = (vi)fa 1 is a bounded solution of the orientable symmetric system (1.9). 
Assume also that either n = 2, 0 < Sj < 1 and v is stable or n = 3, 1/2 < s* < 1 and v is H-monotone. 
Then, there exist a constant T,; £ S'" -1 and v* : R + x R + —> R such that 

(4.10) Vi(x,y) = v*(Ti-x,y) for (x,y)eR” +1 , 

and i = !,••• ,m. Moreover, for all i,j = !,••• ,m vectors V x Ui(x,0) and V x v 3 (x,0) are parallel and the 


angle between two vectors is arccos 


\djHilv)\ 


Hi(v) 


Proof: First, let n = 2 and v be a stable solution of (1.9). We apply Lemma 3.2 and Theorem 4.1 when 
hij = djHi{v)d n Vid n Vj , / is the identity function and F is constant. For each 1 <k< to, set yfax) := Pr(x) 
in the geometric Poincare inequality (3.15) when 

log R, if \x\ < \/R, 

(4.11) Pr(x):={ 2 los f- 1 ° g ' g ' , if \f~R < \x\ < R, 

0, if Ixl > R. 
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From (3.15) and the fact that \Vp R \ < when y/R < x < R, we conclude 


(4.12) 


I log R | 2 I y 1 2si (\Vvi\ 2 Ai+\V Ti \VxVi\\ 2 )dx 

+ |l°gi?| 2 V / (|9ji?i(w)||V x Vj||V x ^| - djHi{v)W x Vi ■ V x Vj)dx 


c 


y l ~ 2si |V«i 


-<dr. 


/s+\B 


R^VR 


Here, we have used the notion of symmetric systems that is yjd 7 H, ( v )<9,; H y (v) = \diHj(v)\. On the other 
hand, Lemma 3.2 implies that for each index 1 < i < m, 


(4.13) [ ?/ 1-2si |Vi>i| 2 dx < CR 2 . 

J K 

Straightforward calculations show that for each i, we have 

(4.14) [ -Ly^Z^Vifdx 

Jb r\ B % 1*1 

= 2 f f T - 3 y 1 - 2s *\S7v i \ 2 dTdx + ^ [ y'-^lVvtfdx 
J B + r \B%J |x| R J B + r \B% 

< 2 [ t~ 3 [ y l ~ 2si \X/Vi\ 2 dxdT + ^— f y 1_2si \X/Vi\ 2 dx. 

Jvr Jb+ R z Jb+ 

Combining this and (4.13), we conclude 

(4.15) [ yl 2Sl }^ Vi][2 dx<ClogR, 

Jb + r \b% M 2 

where the positive constant C is independent from R. The above decay estimates (4.12) and (4.15) imply that 
|Vu i | 2 „4 2 + |Vr i |V, r v i || 2 and \djH i (v)\\V x v i \\V x Vj\-djH i (v)V x Vi-'V x Vj vanish on R+ and 5R+, respectively. 
Therefore, each Vi is a one-dinrensional function and vectors \7 x Vi(x, 0) and V x Vj(x, 0) are parallel and they 
are in the same direction when djHi(v) is positive and in opposite directions when djHi(y) is negative. This 
proves the desired result in two dimensions. 

We now suppose that n = 3 and v is a id-monotone solution of (1.9). Note that id-monotonicity implies 
stability. The fact that v = (vi)tL 1 is a bounded stable solution of (1.9) in R% implies that the function 
v = (vi)il 1 where Vi(xi,X 2 ,y) '■= Hm X3 ^ yoo Vi(xi,X 2 ,X 3 ,y) is also a bounded stable solution for (1.9) in 
R(j_. From our previous arguments regarding R(!_, we conclude reduction of dimension for each u) that is 
Vi(x,y) = h*(r,; • x, y) for (x,y) € R 3 and for some Tj G S 1 . From this and Lennna 2.3, we conclude that 
the energy of v in C R C R+ is bounded by 

(4.16) E r (v) < CR 2 , 

when E r {v) = YrL i 2 T~ Ic R y 12Si Y^i\ 2 dx — J Br {H(v) — eg) dx for eg := sup H(v). Set <Ji ^ when 
</>i := d n Vi and ^ := Vu* • r\ for any fixed y = (?/, 0) G R"^ 1 x {0}. Note that cr = (ai) 7 ^L 1 satisfies (4.6) and 
(4.8). To apply Theorem 4.2, we only need to prove the following energy estimate 


m 1 C 

(4-17) ^2d~ / y 0i l Vv *| 2d * - CR2 X{s,>\/ 2 } +CR 2 \ogR X{ s ,=i/ 2 }- 

i— i Ujs i JCr 

Define the sequence of functions u 4 = (u|)’(L 1 when v\ (x,y) := Vi(x',x n +t,y) for t G R and (x,y) = 
(x',x n ,y) G R+. Note that u 4 is a bounded solution of (1.9), i.e., 

f div(y 1_2si Vu|) = 0 in R” +1 , 

\ — liniy-^o y 1 ~ 2si d y vj = d Si ddj(u 4 ) in <9R" +1 . 


(4.18) 
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Straightforward calculations show that v\ converges to u, in C/ oc (K") for alH = 1, • • • , m and 
(4.19) lim E^v 1 ) — E R (v). 


On the other hand, for every T > 0 we have 

(4.20) E r (v) = E r {v t ) - [ d t E R {y t )dt. 

Jo 

Note that applying similar arguments, as the ones given in the proof of Theorem 3.3, one can get a lower 
bound for dtE^v*) of the form (3.35). From this and (4.20), we have 


(4.21) E r (v) < E r {v t ) + M 


R m l-2s» 

+ M 


ldB R Jo 1 + 2/ 
Now, taking a limit when T —> oo and using (4.16) we get 


I y?/ 2si dx. 

B R y~{y—R} i=1 


(4.22) 


E r (v) < E r (v) + M 


n r-R m i_o o. m 

dB R Jo i=1 1 + 2/ i=1 


,■ ,-R m 1-2Si m 

/ / V ^- dydtW 1 - 1 +MV i?" 

/as«do i=1 ! + y i=1 


< Mi ? 2 + M 

< MR 2 + My J?" 2 s ‘X{o<Si<i/ 2 } + M y i? n 1 X{i/ 2 < Si <i} 

i=l 

m 

+MR n ~ 1 log i? x {si= i/ 2} + M y i?" 


Z =1 
on—2si 


i—1 


Here, we have used estimates provided in Lemma 2.3. Note that when 1/2 < s* < 1 and n = 3, for any 
1 < i < m, we have 


(4.23) n — 2s,; < n — 2s* < n — 1. 

Applying (4.22) when n = 3 and 1/2 < s* < 1, we conclude 

(4.24) E r {v) < MR 2 + Mi ?" -1 < Ci? 2 , 

when C is a positive constant that is independent from R. Similarly, for the case of s* = 1/2 and n = 3, we 
have 

(4.25) E r (v) < MR 2 + Mi ?"” 1 log i? < Ci ? 2 log R. 

This completes the the proof. 

□ 

Methods and Ideas applied in the above proof are strongly motivated by the ones given by Ambrosio and 
Cabre in [2], by Alberti, Ambrosio and Cabre in [1], by Farina, Scuinzi and Valdinoci in [28], by Ghoussoub 
and the author in [31] and Sire and the author in [32]. We now provide another consequence of Theorem 4.1. 
The following theorem clarifies the behaviour of derivatives of each u, in various directions when v = (vi)YLi 
is a bounded stable solution of symmetric system (1.9) in lower dimensions. 


Theorem 4.3. Suppose that v = (vi)YYi a bounded stable solution of symmetric system (1.9) for n < 2 
when s* belong to [ 5 ,1) and for n < 1 + 2s* when s* belongs to (0, 5 ). Then, either each d x Vi(x, y) vanishes 

in R " +1 or it does not change sign in R " +1 for every 1 < i < m. 

Proof: Let v be a bounded stable solution of (1.9). Then, there exits a sequence of functions <f> = (<(>,), such 
that each ^ does not change sign and it satisfies 

J di v{y ai S7f>i) = 0 in R" +1 , 

\ — lim.y^o y ai d y 4>i = d Si J2]LidjHi(v)(t>j in dR" +1 . 
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(4.26) 



For each index 1 < i < to, define the quotient cq := that implies ( (Jifa ) 2 = ( d x Vi ) 2 . From Lemma 3.2, 
we have 


(4.27) 


/ y ai {<Ji4>i) 2 dx < / y ai \\7vi\ 2 dx < CR n 
'Cr Jc r 


1 + 2 / 


dy, 


where i? > 1. Therefore, straightforward calculations show that 


(4.28) 


ICr 



for 

R n In R 

for 

R n 

for 


1 

2 ’ 


2 ’ 


Note that for cq and fa equations (4.6) and (4.8) hold. For n < 2 when all s,; belong to [|,1) and for 
n < 1 + 2s* when at least one of s, belongs to (0, |), estimate (4.3) holds for an appropriate F £ T. Set 
hij := djHi(v)fa(j)j and / to be the identity function in Theorem 4.1. Note that for symmetric systems, 
we have /q ;J = hjj. Theorem 4.1 implies that each cr, is constant. Therefore, there exists a sequence of 
constants C = ((7j)j such that d x Vi(x,y) = Cifa(x,y) for (, x,y ) £ R" +1 . Since each fa does not change sign, 
the proof is completed. 

□ 


The De Giorgi’s conjecture provides a reduction of dimensions, to one-dimension, for bounded monotone 
solutions of the Lane-Enrden equation when n < 8 . The latter theorem provides a counterpart of the 
conjecture to multi-component fractional symmetric systems with a general nonlinearity. In what follows, 
we assume certain extra assumptions on the sign of the nonlinearity H and we establish a Liouville theorem 
for bounded stable solutions of (1.9) in lower dimensions applying Theorem 4.3, Theorem 4.1 and Theorem 
3.1. 


Theorem 4.4. Suppose that v = is a bounded stable solution of (1.9) when either Hifv) > 0 for all 

1 < i < m or Y(!iLi v iHi(v) < 0. Then, each i\ must be constant provided n < 2(1 + s*). 


Proof: Since v is a stable solution, there exists a sequence 4> = {fa) r (L 1 satisfying (1.10). On the other hand, 
applying Lemma 3.3, for each i, we have 


(4.29) f y 1 ~ 2si \Vv i \ 2 dx<CR n - 2si . 

J Cr 

The fact that n < 2 + 2s*, implies that n — 2s,; < n — 2s* < 2. Note that cr,; = v ^' l? satisfies conditions 
of Theorem 4.1 for F(r) = 1, hij = djHi(v)fa(f>j and / to be the identity function. Therefore, each eq 
must be constant for an arbitrary direction r}. This implies that there exist a constant T* € S'™ -1 and 
v* : R + x R + R such that Vi(x,y) = n*(r, ; • x,y) for (x,y) € R" +1 . In other words, v = (n i )+ 1 is a 
bounded stable solution of (1.9) for n = 1. Applying Theorem 4.3, we conclude that d x v, does not change 
sign in R^_. Here, we have used the fact that when d x Vi vanishes, boundedness implies that each must 
be constant. Therefore, iq has to be strictly monotone in x which together with boundedness of proves 
the existence of lim x ^.±oo Vi(x, 0). Let lirn^-oo Vi(x, 0) = li and lim^-nx, Vi(x, 0) = Li, where U and Li 
are constants. Since i>,: is strictly monotone, we conclude L < Li. We now apply the Hamiltonian identity 
provided in Theorem 3.1 when x —> ±oo to get H(l) = H(L ) where l = {li) r )L 1 and L = {Li) r fL 1 . Note that 
this is in contradiction with the following 


(4.30) 


m 

0 = H{L) - H(l) = ^( Li - k)Hi(t(L -l) + /), 


for some t £ (0,1). This completes the proof. 

□ 

Mathematical techniques and ideas that we applied in the above proof are strongly motivated by the ones 
given in [23] by Dupaigne and Farina. In [23], authors proved that any bounded stable solution of (1.12), that 
is when m = 1 and s = 1 in (1.1), is constant provided n < 4 and 0 < H £ C' 1 (R) is a general nonlinearity. 
Note that for particular nonlinearities the critical dimension is much higher than four dimensions. We also 
refer interested readers to [5] by Cabre and Capella, to [49] by Villegas and to [30] by the author for the case 
of radial stable solutions where the optimal dimension is n = 10 for a general nonlinearity H £ C 1 (R). So, we 
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expect that Theorem 4.4 could be improved. For specific nonlinearities H{u) = e“, H(u) = u p where p > 1 
and H(u) = —u~ p where p > 0 the equation is called Gelfand, Lane-Emden and Lane-Emden with negative 
exponent equations, and Lioville theorems are given for the following optimal dimensions, respectively, 

• 1 < n < 10 by Farina in [27], 

• 1 < n < 2 + + \/p[p — 1)) by Farina in [26], 

• 1 < n < 2 + y/p(p + 1)) by Esposito, Ghoussoub and Guo in [25]. 

In the next section, we study a two-component nonlinear Sclrrodinger system, that is a particular case of 
(1.1), and we prove various Liouville theorems. 


5. Nonlinear Schrodinger System; Particular Nonlinearity 


Consider the following two-component system 
, , f (—A) Sl ui = piu\ + f3u\ui in R n , 

\ (—A Y 2 U 2 = P 2 ul+Pulu 2 in R”, 

when 0 < si,S 2 < 1 and AWATj/? are parameters. This is a special case of system (1.1) when m = 2 and 
Hi(u\,U 2 ) = ATuf + (3u\u\ and H 2 (u\,U 2 ) = AT+ / 3u\u2■ The above system arises in Bose-Einstein 
condensations and it is well-studied in the literature. We refer interested readers to [3, 15, 32, 43, 46, 51-54] 
and references therein for more information. The extension function pair (vi,V 2 ) given in (1.9) satisfies 


(5.2) 


f di v(y ai X7vi) = 0 in R” +1 , 

\ —d Si lim^o y ai d v Vi = H i (v 1 ,v 2 ) in 9R” +1 , 


Q2s i — 1 p / \ 

when a, = 1 — 2s* and d Si = - f° r 1 < i < 2. Note that when (3 = 0, system (5.1) becomes decoupled 

and each equation in (5.1) is of the from 

(5.3) (-A ) a w = w 3 in R n . 


when /ii,/i 2 > 0 for either w = yjJT[ui and a = s i or w = yfp 2 U 2 and a = S 2 - The above equation (5.3) 
is known as the fractional Lane-Emden equation and nonnegative solutions of this equation are classified 
completely in the literature. It is known that whenever 


(5.4) n < 4a, 

the only nonnegative solution for (5.3) is the trivial solution and n = 4a is the critical dimension, see 
[14, 38, 42]. We refer interested readers to [16, 17, 33] for the classification of stable solutions of (5.3) when 
0 < a < 2. In this article, we are interested in the case of /3 ^ 0. The following Liouville theorem addresses 
stable solutions of (5.2) and is a direct consequence of Theorem 4.4. The proof is straightforward and we 
omit it here. 


Theorem 5.1. Let v = (v\,V 2 ) be a bounded stable solution of (5.2) when n < 2 + 2min{si,S2j-- Assume 
that either AT, ^2 < 0 and |/3| < y/JLifLf or at, a l 2iP > 0 and solutions (iq, V 2 ) are nonnegative. Then, each 
Vi = Ci where Ci is constant. 

We now provide a Liouville theorem for solutions of (5.1) in the absence of stability. 

Theorem 5.2. Suppose that s\ = S 2 = a and AT,AT are nonnegative. Assume that u = (u\,U 2 ) is a 
nonnegative solution of (5.1) when n < 3a and f3 > then (u\,U 2 ) = (0,0). 

Proof: Let u = (u\,U 2 ) be a nonnegative solution of (5.1) when si = S 2 = a. Without loss of generality 
assume that u\ > 0 and at > 0. It is straightforward to show that equivalently {u\,U 2 ) is a solution for the 
integral system 


f Mx) = f Rn ^^l Ul{v) dy former, 

\ u 2 (x) = f Rn M2U2 (v) dy for x G R n 

Assume that parameter (3 is nonnegative. Since each iq is nonnegative, system (5.5) gives 

f U!(x) > at f R n j x l‘yin-^ dy for x G R", 

l U ^ X ) ^ 


(5.5) 


(5.6) 


“■ 2 ,I ( n-2» dy for x G 
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Applying Liouville theorems given in [12], see also [13, 14, 42], to each of the above integral inequalities we 
conclude that (ui, u 2 ) = (0,0) when n < 3a. We now let /3 be negative that is —y/fj, 1 H 2 < /? < 0. Multiplying 
the first equation in (5.5) with fj p, 2 / and adding both equations in (5.5) we get 

(5.7) 


\Jfi 2 /Hiu\{x) + u 2 (x) = 

f 1 


^W/h (mu\{y) + fiu!(y)ul(y)) + p 2 u\{y) + Pu 2 {y)u\(y) 


dy , 


for x € 

M such that 


J K» \x - y\ n ~ 2a L 

. To simplify the right-hand side of above equality, we claim that there exists a positive constant 


(5.8) v p 2 / Hi (/-ti u\ + Pu\u 2 ) + p 2 u\ + fiu 2 u\ > M y \/p 2 /piUi + u 2 
In order to prove this claim, let r := u 2 /u\ and define continuous function F : R + — 

(5.9) F{r) = 


as 


\//i2/Vi (/h + (3t 2 ) + p 2 r 3 + /3t 


(v^W/h +t ) 3 

Note that F( 0) = > 0 and lim,--^ F(t) = p 2 > 0. In addition, for (3 > — ^/pip 2l we have 

(5.10) (V/ 12 /M 1 + r ) F(t) > \/V 2 /Vi (ah - v'MiM 2 T 2 ) + t(^ 2 t 2 - 

(5.11) = (v 7 /^ 7- - v'Tl) 2 + \Zhih2) > 0, 


that is F(t) > 0 when t > 0. This proves the claim. From (5.8) and (5.7), we get 
---- /' 1 / -- \ 3 

(5.12) A/ ~ ' ' ' ' ' ” 


\J 112 /HiUiix) + u 2 {x) > M 


In— 2a. 


/\x - y\ 

Let z(x) := \[M p 2 /piu\(x) + u 2 (x)^ > 0 in (5.12) to obtain 

z{x) > [ 

Jr 


(v^/Ah ui{y) + u 2 (y)^ dy for 


X G 


(5.13) 


z 3 (y ) 


/r» k - y | ? 


for x G 


We now apply Liouville theorems given in [12, 13] for the above integral inequality to conclude that z = 0 
when n < 2a and for n > 2a when 3 < ff 2a . This implies that whenever n < 3a, we have z = 0. This 
completes the proof. 

□ 

According to (5.4), one may expect that the critical dimension for the above theorem is n = 4a and this 
remains an open problem. Note that methods and ideas applied in the above proof are strongly motivated by 
the ones given in [15] by Dancer, Wei and Wetlr and in [43] by Lin and Wei. We now provide a monotonicity 
formula for solutions of the two-component Schrodinger system (5.2). This is a direct consequence of Theorem 
3.4 and we omit the proof. Note that Frank and Lenzmann in [34] and with Silvestre in [35] used similar 
monotonicity formulae to study uniqueness of solutions for the fractional Schrodinger operator. 

Theorem 5.3. Let v = (vi,v 2 ) be a bounded solution of (5.2). Suppose that pi,p 2 < 0 and |/3| < y/pip 2 . 
Then 


(5.14) 


I{R) = R 


—n+2s 


Y. ds 


y ai \Wvi\ 2 dx - 


J B R x{y= 0} 


(y « 1 + y *4 + pv\vi) dx 


is a nondecreasing function of R> 1. 

Lastly, we provide a monotonicity formula for radial solutions of (5.2) for all parameters /ii,/i 2 ,/3 £ R. 

Theorem 5.4. Let v = (vi)'fL 1 be a bounded radial solution of (5.2), i.e. Vi(x,y) = Vi(\x\,y). Then , the 
following function is nondecreasing in r > 0, 


(5.15) 


J {r) ■■= J2d s 


y 1 2si [(d r Vi{r,y)) 2 - (d y Vi(r,y)) 2 ] dy - Ff(ui(r, 0), u 2 (r, 0)), 
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where H(y) = ^vf + fAv f + More precisely, 


(5.16) 


1 2 «oo 

r 'H Jo 


Proof: Suppose that i> = (i>j) is radially symmetric in x. For r = \x\ we have 

drrVi + n ^d r v i + d yy Vi + Y<Jy v i = 0 in (0, 00) x (0, oo), 

-d Si lim y ^o V ai d y Vi = Hi{v) in (0, oo) x (y = 0), 
where ai = 1 — 2si. Define the following function of r, 

2 7 /*no 


(5.17) 


(5.18) 


P(r) ■= Y ' 


■y 1 2si [(.d r Vif - {dyVi) 2 ] dy 


i=l 


Taking derivative of p with respect to r and using (5.17) to substitute values of d rr Vi , we conclude 

_ i 2 «oo 2 /-oo 

(5.19) p'(r) = -- y ai (d r Vi) 2 dy ~y2d Si y ai (d y yVi)(d r Vi)dy 

r i=1 J o i=1 do 

2 «oo 2 /* oo 

-Y d »t a i y~ 2si {drVi)(d y v z )dy - Y ds * / y ai ( d y v i)(dr y Vi)dy. 

i= 1 do i = l do 

Applying integration by parts and using the boundary term in (1.9), we have 

2 /»oo 2 

(5.20) -V'ds, / [y ai <9 TO 't>j<9 r Vj + cny~ 2si d r VidyVi + dy = lim V] d Si y ai d r Vid y Vi 

■ i do , 

%— 1 2=1 

= d r (H(v)Y 

Combining (5.20) and (5.19) and setting J(r) := p(r) — H(v(r, 0)) completes the proof. 


□ 


We end this section with this point that in the absence of monotonicity and stability assumptions, the 
qualitative behaviour of solutions of elliptic and Hamiltonian systems with a general nonlinearity are studies 
extensively in the literature. We refer interested readers to [18, 19, 39] and references therein. 
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